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Abstract 

This paper is the first in a series revisiting the Faraday effect, or more 
generally, the theory of electronic quantum transport/optical response in 
bulk media in the presence of a constant magnetic field. The independent 
electron approximation is assumed. At zero temperature and zero fre- 
quency, if the Fermi energy lies in a spectral gap, we rigorously prove the 
Widom-Streda formula. For free electrons, the transverse conductivity 
can be explicitly computed and coincides with the classical result. In the 
general case, using magnetic perturbation theory, the conductivity ten- 
sor is expanded in powers of the strength of the magnetic field B. Then 
the linear term in B of this expansion is written down in terms of the 
zero magnetic field Green function and the zero field current operator. 
In the periodic case, the linear term in B of the conductivity tensor is 
expressed in terms of zero magnetic field Bloch functions and energies. 
No derivatives with respect to the quasi-momentum appear and thereby 
all ambiguities are removed, in contrast to earlier work. 



1 Introduction 

In sharp contrast with the zero magnetic field case, the analysis of properties of 
electrons in periodic or random potentials subjected to external magnetic fields 
is a very challenging problem. The difficulty is rooted in the singular nature 
of the magnetic interaction: due to a linear increase of the magnetic vector 
potential, the naive perturbation theory breaks down even at arbitrarily small 
fields. 

To our best knowledge, only the periodic case has been considered in con- 
nection with the Faraday effect for bulk systems. The first full scale quantum 
computation was done by Laura M. Roth jSJ (for a review of earlier attempts 
we direct the reader to this paper). The physical experiment starts by sending 
a monochromatic light wave, parallel to the Oz direction and linearly polarized 
in the plane xOz. When the light enters the material, the polarization plane can 
change; in fact, there exists a linear relation between the angle 9 of rotation of 
the plane of polarization per unit length and the transverse component of the 
conductivity tensor axy (see formula (1) in The material is chosen in such 
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a way that when the magnetic field is zero, this transverse component vanishes. 
When the magnetic field B is turned on, the transverse component is no longer 
zero. For weak fields one expands the conductivity tensor to first order in B 
and obtains a formula for the Verdet constant. 

Therefore the central object is axy{B), which depends among other things on 
temperature, density of the material, and frequency of light. Using a modified 
Bloch representation. Roth was able to obtain a formula for ^^^^ (0), and studied 
how this first order term behaves as a function of frequency, both for metals 
and semiconductors. 

However, the theory in [2] is not free of difficulties. First, it seems almost 
hopeless to estimate errors or to push the computation to higher orders in B. 
Second, even the first order formula contains terms which are singular at the 
crossings of the Bloch bands. Accordingly, at the practical level this theory 
only met a moderate success and alternative formalisms have been used, as for 
example the celebrated Kohn-Luttinger effective many band Hamiltonian (see 
^lEHlEHl and references therein), or tight-binding models Since all these 
methods have limited applicability, a more flexible approach was still needed. 

In the zero magnetic field case, a very successful formalism (see e.g. |27II16I 
12 E] and references therein) is to use the Green function method. This is based 
on the fact that the traces involved in computing various physical quantities 
can be written as integrals involving Green functions. The main aim of our 
paper is to develop a Green function approach to the Faraday effect, i.e. for the 
conductivity tensor when a magnetic field is present. Let us point out that the 
use of Green functions (albeit different from the ones used below) goes back at 
least to Sondheimer and Wilson 33 in their theory of diamagnetism of Bloch 
electrons. Aside from the fact that the Green function (i.e. the integral kernel 
of the resolvent or the semi-group) is easier to compute and control, the main 
point is that by factorizing out the so called "non integrable phase factor" (or 
" magnetic holonomy" ) from the Green function, one can cope with the singular- 
ities introduced by the increase at infinity of the magnetic vector potentials. In 
addition, (as it has already been observed by Schwinger 1221 in a QED context) 
after factorizing out the magnetic holonomy one remains with a gauge invariant 
quantity which makes the problem of gauge fixing irrelevant. The observation 
(going back at least to Peierls [2^1) that one can use these magnetic phases 
in order to control the singularity of the magnetic perturbation has been used 
many times in various contexts (see e.g. (HSIEI)- We highlight here the re- 
sults of Nedoluha [21]] where a Green function approach for the magneto-optical 
phenomena at zero temperature and with the Fermi level in a gap has been 
investigated. 

But the power of this method has only recently been fully exploited in 
jl3l lllj , and developed as a general gauge invariant magnetic perturbation the- 
ory in |23j . Applied to the case at hand, this theory gives an expansion of the 
conductivity tensor in terms of the zero magnetic field Green functions. More- 
over, it is free of any divergences. A key ingredient in controlling divergences is 
the exponential decrease of the Green functions with the distance between the 
arguments, for energies outside the spectrum ^||221- We stress the fact that 
since no basis is involved, periodicity is not needed and the theory can also be 
applied to random systems. Finite systems and/or special geometries (layers) 
are also allowed. The content of the paper is as follows: 

In Section 121 we give a derivation of the conductivity tensor from first prin- 
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ciples in the linear response theory. We include it to point out that it coincides 
with various formulas used before. Although in physics establishing the Kubo 
formula is considered somehow a triviality, from a mathematical point of view 
it remains a serious chalenge (see |H])- 

Section|3contains the precise formulation of the thermodynamic limit, stated 
in Theorem 13. II We do not give its proof here, but we try to explain why it is 
true. 

Section^shows that at the thermodynamic limit, at zero temperature, zero 
frequency, and for the Fermi energy in a spectral gap, we re-obtain a formula of 
Streda [31] for the transverse component of the conductivity tensor, known from 
the Integral Quantum Hall Effect (IQHE). A precise statement and its proof 
are contained in Theorem 14. II Moreover, under the proviso that exponentially 
localized Wannier function exist (see Theorem 14. 2|l . this transverse component 
vanishes (see also [HHl; |H] for related results). We stress that this result holds 
for the whole axy{B) as long as the magnetic field is not too large, not just 
for ^^^(0). The vanishing of its first order correction was in fact claimed in 
formula (50) in j31!. 

Section contains the exact quantum computation of axy (B) for free elec- 
trons; in spite of the fact that such a result might be known (and it is known 
at zero frequency), we were not able to find it in the literature. Interestingly 
enough, the quantum computation gives the same result as the well known 
classical computation (when the relaxation time is infinite). 

Section El contains the core of the paper, which includes the derivation of 
^^^(0) for general Bloch electrons. As in the zero magnetic field case, its 
formula only contains zero magnetic field Green functions and current operators. 

Section 13 deals with periodic systems, and the result of the previous section 
is written down in terms of zero magnetic field Bloch functions and bands. 

At the end we have some conclusions. 

The main goal of this paper is to present the strategy, state the results 
concerning the Verdet constant, and to outline future theoretical and practical 
problems. Detailed proofs of the thermodynamic limit and of other technical 
estimates will be given elsewhere. 

2 Preliminaries. The conductivity tensor in the 
hnear response regime 

We begin by fixing the notation used in the description of independent electrons 
subjected to a constant magnetic field. The units are chosen so that h = 1. Since 
we consider spin 1/2 particles, the one particle Hilbert space for a non-confined 
particle is 

Woo =i^(K^)©i2(R3) 

with the standard scalar product. Accordingly, all operators below and their 
integral kernels are 2x2 matrices in the spin variable. We choose the constant 
magnetic field of strength B to be oriented along the z-axis. Then the one 
particle Hamiltonian with the spin-orbit coupling included is (see e.g. |31|') 

H^{B)^^PiB)^ + V + g^IbBa3, (2.1) 
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with 

P(B) = -iV -ba+ -^s A (VV) = P(0) - 6a (2.2) 

where 

b=-'-B 
c 

and a(x) is an arbitrary smooth magnetic vector potential which generates a 
magnetic field of intensity B = 1 i.e. V A a(x) = (0, 0, 1). The most frequently 
used magnetic vector potential is the symmetric gauge: 

ao(x) = in3Ax. (2.3) 

where 113 is the unit vector along z axis. 

In the periodic case we denote by C the underlying Bravais lattice, by Q its 
elementary cell and by O* the corresponding Brillouin zone. and stand 
for the volumes of the elementary cell and Brillouin zone respectively. In the 
absence of the magnetic field one has the well known Bloch representation in 
terms of Bloch functions: 

*,(x,k) = -^=e'^''u,(x,k), X G (2.4) 

where Uj (x, k) are the normalized to one eigenfunctions of the operator 

h{\i)uj (x, k) = \j {k)uj (x, k) (2.5) 

'^('^)=2^(-^^^+2i?^^(^^) + '^)' + ^' 

= ^(p + k)^ + V, kef]*, (2.6) 

P = -^^-+2i?^^(^^)' 

defined in L^(0) © L^(0) with periodic boundary conditions. We label Aj (k) 
in increasing order. We have to remember that, as functions of k, Aj(k) and 
Mj(x, k) are not differentiable at the crossing points. Since the \I'j(x, k)'s form 
a basis of generalized eigenfunctions, the Green function (i.e. the integral kernel 
of the resolvent) writes as: 

GS'(x,y;.)=/ E^^^lf^*. (") 

and it is seen as a matrix in the spin variables. The above formula has to be 
understood in the formal sense since the series in the right hand side is typically 
not absolutely convergent, and care is to be taken when interchanging the sum 
with the integral. Notice however that Gtc^(x, x';2:) is a well behaved matrix 
valued function. 

We consider a system of noninteracting electrons in the grand-canonical 

ensemble. More precisely, we consider a box Ai C M^, which contains the 
origin, and a family of scaled boxes 

Az, = {x e : x/L e Ai}. (2.8) 
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The thermodynamic hmit will mean L oo, that is when tends to fill out the 
whole space. The one particle Hilbert space is := L^(Al)0L^(Al). The one 
particle Hamiltonian is denoted by Hl{B) and is given by 1)2. l|l with Dirichlet 
boundary conditions (i.e. the wave- functions in the domain of Hl{B) vanish at 
the surface OAl). More precisely, we first define it on C^{Al) © C^{Al), and 
then Hl{B) will be the Friedrichs extension of this minimal operator. This is 
indeed possible, because our operator can be written as (up to some irrelevant 
constants) —Auh + W, where Ad is the Dirichlet Laplacian and is a first 
order differential operator, relatively bounded to —Ad/2 (remember that L < 
oo) with relative bound zero. The form domain of Hl{B) is the Sobolev space 
Hq{Al) © Hq{Al), while the operator domain is 

Dom(i/L(S)) =Dl(B Dl, Dl := H^{Al) n HI{Al). (2.9) 
Moreover, Hl{B) is essentially self-adjoint on Cj^^(Al), i.e. functions with 

support in A^ and indefinitely differentiable in A^, up to the boundary. 

We assume that the temperature T = l/(fc/3) and the chemical potential /i 
are fixed by a reservoir of energy and particles. We work in a second quantized 
setting with an antisymmetric Fock space denoted by J^l ■ Denote the operators 
in the Fock space with a hat and borrow some notation from the book of Bratelli 
and Robinson [Jj: if A is an operator defined in TYl, we denote hy A = dT{A) 
its second quantization in the Fock space. At t = —qg the system is supposed 
to be in the grand-canonical equilibrium state of temperature T and chemical 
potential /x, i.e. the density matrix is 

^ e-'^*", (2.10) 



Tr(e-/3^-) 
where 

k^=dT{HL{B)-ii-ld) (2.11) 

is the "grand-canonical Hamiltonian" . 

The interaction with a classical electromagnetic field is described by a time 
dependent electric potential 

F(x, t) := (e*"* + e-^^*)eE • x, i < 0, x e Al. (2.12) 

so the total time dependent one-particle Hamiltonian is 

H{t)^HL{B) + V{t). (2.13) 

Notice that e near E is the positive elementary charge. Here we take Im a; < 
which plays the role of an adiabatic parameter, and insures that there is no 
interaction in the remote past. Finally, the one-particle current operator is as 
usual 

J = -ez[i/i(B),X] ---P(S), (2.14) 
m 

where X is the multiplication by x. Note that J is a well defined operator on the 
domain of Hl{B), because multiplication by any component of X leaves this 
domain invariant (see H2.9(l ^. Moreover, since L < cxd, X is a bounded operator. 
In fact, X is the true physical self-adjoint observable, while P(-B) (or J) appear 
when one differentiates the map 1 1-^ gitHL{B)-^^-itHL(B) ^j^g strong sense on 
the domain of Hl{B). 
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We assume that the state of our system is now described by a time-dependent 
density matrix, p{t), obtained by evolving from — oo up to the given time, 
i.e. 

idtp[t) = [H{t),p{t)], p(-oo) = pe. (2.15) 

Going to the interaction picture and using the Dyson expansion up to the first 
order, one gets 

p{t = Q) = pe-i [dT{V{s),pe]ds + 0{E^), (2.16) 

J —oo 

where 

V(s) ■.^e"'"^^^^V{s)e-''^^^^l (2.17) 
The current density flowing through our system at i = is given by (see (|2.16|l '): 

j = ^Tr,, (p(O)j) ^ J^Tr,, (^.j) 

' Tr^, [dT{V{s)),p,]3ds^ +0(E^). (2.18) 



|A 

In evaluating the r.h.s. of (|2.18() we use the well known fact that traces over the 
Fock space can be computed in the one-particle space (see Proposition 5.2.23 in 

my- 

Tr^, {ped^iA}} = Tin, {fFD{HL{B))A} , (2.19) 
where fpD is the Fermi-Dirac one-particle distribution function: 

fpoix) := xeR,P>0,peR. (2.20) 

Plugging ifm^ into ^J^, the identity [dr(A), dr{B)] ^ dr{[A, B]), the invari- 
ance of trace under cyclic permutations and ignoring the quadratic correction 
in E one arrives at 

i = r^^TTnAfFDiHL{B))J} 

\^l\ 



I e 
|Al| to 



Trw. [V{s),'P{B)]fFn{HL{B))ds^ . (2.21) 



The first term in (|2.2HI is always zero because of the identity (trace cyclicity 
again) 

Tr„, {[HL{B),X]fFn{HL{B))} = Tr«, {[fFD{HL{B)), HUB)]^} = 0- 

(2.22) 

which is nothing but the fact that the current vanishes on an equilibrium state. 
Note that these operations under the trace sign are quite delicate, since un- 
bounded operators are involved. Let us for once give a complete proof to H2.22|l . 
We have the identity between bounded operators (consider the first component 

[HL{B),X{\fFD{HL{B)) = HL{B)X^fFD{HL{B)) - XiHL{B)fFD{HL{B)). 

(2.23) 
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Remember that Xi is a bounded operator in the box, and preserves the domain 
of Hl{B). This means that the operator Ol = {Hl{B) + i)Xi{HL{B) + is 
bounded. Hence we can write 

HL{,B)X^fFD{HL{B)) = [1 - i{Hl{B) + i)-^] Ol [HUB) + i]fFD{HL{B)). 

Now the operator [Hl{B) + i\fFD{HL{B)) stiU is trace class due to the expo- 
nential decay of fpD, while [1 — i{HL{B) + i)^^] and Ol are bounded. Thus 
HL{B)XifpD{HL{B)) is trace class and we can compute its trace using the 
complete eigenbasis of Hl{B), which gives the same result as for the other 
operator XiHL{B)fFD{HL{B)). Thus I^TTB) is proved. 
Using ()2.12|l and l|2.17|l one can write 

3 

i« = + '^c.i3{-uj)}Efi, a e {1, 2, 3}, < 0, (2.24) 

0=1 

where the conductivity tensor is given by 

aa.p{B,u)^ (2.25) 
-Tr«, / [e''"-^''^xpe'''''-'^''\PUB)]fpD{HL{B)y'''ds. 



i e2 



I Ail m 



Performing an integration by parts, using the formulas i[HL{B), xp] = PpiB) /m 
and i[Pa{B), x/j] = Sap one arrives at 

a^piB,Lu) = -^^{S^^,TiifpDiHLiB))) (2.26) 
\Al\ imuj 

+ -TrT e^^("+«-(^))p^(i?)e-^^^-(^)[P„(i?),/fD(i?L(i3))]ds}, 
■m- J-oo 

and this coincides (at least at the formal level) with formula (5) in [SJ. Notice 
that from now on, we write just Tr when we perform the trace, since we only 
work in the one-particle space. 

Since we are interested in the Faraday effect, and we assume that the mag- 
netic field B is parallel with the z axis, we will only consider the transverse 
conductivity ai2{B,uj). Hence the first term vanishes. We now perform the 
integral over s with the help of Stone's formula followed by a deformation of the 
contour (paying attention not to hit the singularities of fpL)(z) or to make the 
integral over s divergent 

fpDiHLiB))e'^^"-^^^+^^ = / fpD{z)e''^''+^\HLiB) ~ z^dz. (2.27) 

where 77 is either or w, the contour is counter-clockwise oriented and given by 
T^^{x±id: a <x < oo]\J {a + iy: -d < y < d} (2.28) 

with 

TT llm Ll>\] _ 
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and a + 1 lies below the spectrum of Hl{B). As a final result one gets 



ai2(i3,^) = -- ^-r- (2.30) 

• Tr j fpoiz) {Pi{B){Hl{,B) - z)-^ P2{B){Hl{B) - z - uo)-^ 

+ z— >z — cj}dz=: — —a^fB, uj) 

where "z — > z — cj" means a similar term where we exchange z with z — w. 
Now one can see that by inserting the eigenbasis of Hl{B) one obtains the well 
known formula derived from semi-classical radiation theory (see e.g. formula 
(4) in ^). 

3 Gauge invar iance and existence of the ther- 
modynamic Umit 

Up to now the system was confined in a box A/,. As is well known (see e.g. 
a direct evaluation of H2.3()|l (or previous formulas equivalent to it including 
formula (4) in Roth's paper) is out of reach: the eigenvalues and eigenstates 
of H{B) are rather complicated (even in the thermodynamic limit A^ — > K^) 
and at the same time the Bloch representation is plagued by singular matrix 
elements of the magnetic vector potential. Roth used a modified magnetic Bloch 
representation in jHOI and derived a formula for the linear term in B of H2.3()|l in 
terms of the zero magnetic field Bloch representation. Still, her procedure is not 
free of difficulties since it involves VkUj(x, k) which might not exist at crossing 
points. In addition, it seems almost hopeless to control the errors or to push 
computations to the second order in B which would describe the Cotton-Mouton 
effect for example. 

In what follows, we shall outline another route of evaluating H2.30|l which is 
mathematically correct, systematic, and completely free of the above difhcul- 
ties. There are two basic ideas involved. The first one (going back at least to 
Sondheimer and Wilson in their theory of diamagnetism) consists in writing 
the trace in H2.30|l as integrals over A^ of corresponding integral kernels. This 
is nothing but the well known Green function approach (see e.g. |15| ) which 
has been very successful in computing optical and magneto-optical properties of 
sohds (see e.g. JHl, |2j) in the absence of an external magnetic field. The 
point is that the integral kernels are on one hand easier to control and compute, 
and on the other hand they do not require periodicity. Moreover, this approach 
proved to be essential in deriving rigorous results concerning the diamagnetism 
of free electrons ^] and actually we expect the methods of the present paper 
to simplify the theory of diamagnetism of Bloch electrons as well. 

However, when applying Green function approach in the presence of an exter- 
nal magnetic field one hits again the divergences caused by the linear increase of 
the magnetic vector potential: naively, at the first sight ai(i?, uj) is not bounded 
in the thermodynamic limit L —> oo but instead grows like the second power of 
L. It was already observed in that these divergent terms vanish identically 
due to some identities coming from gauge invariance. 
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This is indeed the case and the main point of this paper is to show, following 
the developments in d]i |22|> that factorizing the so called " non-integrable 
phase factor" from the Green function (the integral kernel of {Hl{B) — C)~^) 
allows, at the same time, to eliminate the divergences coming from the increase 
of the magnetic vector potential and to obtain a controlled expansion in powers 
of B. In addition, this leads to expressions of aL[B,u}) which are manifestly 
gauge invariant. 

For an arbitrary pair of points x, y G consider the "magnetic phase" 
associated with the magnetic vector potential a(u) defined as the path integral 
on the line linking y and x: 

'/>a(x,y) = / a(u) • du. (3.1) 

The magnetic phase satisfies the following crucial identity: for every fixed c 

g-^b0.(x.c)p(^^g^60^(x,c) ^ p^Q^ _ ^^^^ _ ^3 3) 

where A(x) = A x, i.e. irrespective of the choice of a(x), 

A(x - c) = ing A (x - c) (3.3) 

is the symmetric (transverse, Poincare) gauge with respect to c. 

Write now the Green function (as a 2 x 2 matrix in the spin space) 

GL(x,y;C) = (i/L-C)"'(x,y) (3.4) 

in the factorized form 

GL(x,y;C) = e'^^»(-'y)ifL(x,y;C). (3.5) 

It is easy to check that while GL(x,y;^) is gauge dependent, i^i(x, y;(^) is 
gauge independent i.e. the whole gauge dependence of Gl(x, y;C) is contained 
in the phase factor e'^'^-'^^^^^'. Plugging the factorization (|3.5|) into the integrand 
of the r.h.s. of (|2.30() . using (|3.2|) and H3.3|l . one obtains that its integral kernel 
writes as 

^t,(x,x')=e'^^-('''''') (3.6) 

2 

■/ dzfFD{z)y^[ dye'^*('''''--'){[(Pi,x(0)-5Ai(x-y))ifi(x,y;z)],,, 

• [(^'2,y(0) - 6A2(x - y))KL{y, x'; z + uj)]^^^' +z^z-uj}, 
where 

$(X, y, x') = 0a(x, y) + 0a (y, x') + <?!)a(x', x) 

is the flux of the magnetic field (0,0, 1) through the triangle A(x, y,x'). Now 
the fact that there are no long range divergences in the formula for ^^^^'(x, x') 
follows from the exponential decay of Green functions jJOl (see also |23): for C 
outside the spectrum of H there exists to(C) > such that as |x — y| ^ oo 

|i^L(x,y;C)| = |GL(x,y;C)| ^e-™(C)l--''l. (3.7) 
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It can be proved (the technical details which are far from being simple will 
be given elsewhere) that (x, x') is jointly continuous and moreover outside 
a thin region near the surface of one can replace it by the integral kernel 
(x, x') of the corresponding operator on the whole M^. Accordingly, up to 
surface corrections: 

1 ^ r 

aL{B,u;)^-——Y, A°°.(x,x)(ix. (3.8) 

Notice that due to the fact that $(x, y,x) = (/)a(x, x) = the phase factors 
appearing in (|3.6|l reduce to unity in H3.8|l . 

In the periodic case, from the fact that in the symmetric gauge the Hamilto- 
nian Hoo{B) commutes with the magnetic translations (actually one can define 
magnetic translations for an arbitrary gauge, just first make the gauge trans- 
formation relating a(x) to A(x)) generated by £, it follows that for 7 G £ we 
have: 

i^oo(x + 7,y + 7; C) = K^{yL,y;C), 

which implies that 

(x + 7, X + 7) = AZ (x, x) 
is periodic with respect to C, hence up to surface corrections: 

1 ^ f 

ai(B,^)«a(B,c^) = --— ^ / ^^,(x,x)dx. (3.9) 

""I I Jn 

Therefore, the transverse conductivity writes as 

ai2iB,Lo)^^a{B,Lu) (3.10) 

with a{B,uj) given by the r.h.s. of H3.9|l . 

A precise formulation of this result is contained in the following theorem: 

Theorem 3.1. Assume for simplicity that is the unit cube in M.^ . The above 
defined transverse component of the conductivity tensor admits the thermody- 
namic limit; more precisely: 

i. The following operator defined by a B{L^ ® L'^)-norm convergent Riemann 
integral, 

Fl ■■= I fFD{z){Pi{B){HL{B) - z)-^P2{B){Hl{B) - z - lu)-' 

+ z^z~Lu}dz, (3.11) 

is in fact trace-class, and crj2^(i?,a;) = m^t2|Ai, | ^^(Fl)- 

ii. Consider the operator Foo defined by the same integral but with Hoo{B) 
instead of Hl{B), and defined on the whole space. Then Foo is an integral 
operator, with a kernel A'^^, (jx., x.' ) jointly continuous on its spatial variables. 

Moreover, the function defined by M.'^ B x ^ sb(x) :— -^s°s(^' ^) G C is 

continuous and periodic with respect to "li^ . 

iii. The thermodynamic limit exists: 

a^^\B,uj) := hm a^^\B,oj) = / SB(x)dx. (3.12) 
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The proof of this theorem wiU be given elsewhere . 



4 The zero frequency limit at T — 0: a rigorous 
proof of the Widom-Streda formula for semi- 
conductors 

Doing some very formal computations one can show that at T = and uj ~ 0, 
(Ti2 {B, Lu) as given by (|3.9(l and (|3.1Q(I coincide with the formula for the quantized 
Hall conductivity (see e.g. formulas (5), (6) in [^) which in turn gives (again 
at the heuristic level) the well known Widom-Streda formula. The original 
derivation has little mathematical rigor, in particular because it assumes some 
very strong assumptions on the existence and regularity of (HooiB) — X+iO)~^ as 
a function of A G K. These assumptions are clearly not true in many situations. 

Here we will show how the Widom-Streda formula can be rigorously obtained 
when the Fermi energy lies in a spectral gap. The problem in which the Fermi 
energy is in the spectrum, remains open. Now assume that for some B, the 
chemical potential /x lies in a spectral gap of Hoo (B). More precisely, throughout 
this section we suppose that {di,d2) C p{Hoo{B)) with di < d2, and take 
/i G {di,d2)- For simplicity, assume that /i = This is the typical situation 

for semiconductors and/or isolators. In the absence of spin, the Widom-Streda 
formula roughly states that: 



dB 



(4.1) 



where N{B,E) is the integrated density of states up to the energy E. When 
the spin is present (this was not considered by Streda), this formula is slightly 
changed. If we denote by Bi the B multiplying the spin matrix ag, in our 
Hamiltonian 1)2.1(1 . and with B2 the B near A, then in fact we have 



dN(Bi,B2,E) 



dB2 



(4.2) 

E=fj., Bi=B2=B 



In the rest of this section we give a rigorous (but still not fully technical) 
proof of (|0) . 

Theorem 4.1. Consider the conductivity at the thermodynamic limit given in 
(I3.12() . and drop the superscript 00. Then if we first take the limit T \ 0, and 
after that — > 0, we get: 



d 1 ^ f 

lim|hnai2(i?,r,u;) = ec— — ^ H^,(x, x)dx 



, (4.3) 

Bi —B2 —B 



where 



Bi I 1 



TV" ^— dz (4.4) 

with a positively oriented contour T enclosing the spectrum of Hoo{B) below fj,, 
i.e. H^ is the Fermi projection onto the subspace of "occupied" states atT = Q. 
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Remarks: 



Streda did not consider spin in his work j34| and in this case the derivative 
with respect to the magnetic field appears in the r.h.s. of (|4.3(l . 

While it is not clear that n^(x, x) is well defined ((iJoc(-B) — z)^^(x, x) 
does not exist!) this can be seen by writing for some a G p{Hoo{B)): 



27r 7r 



27r 



(z - a){Hoo{B) - z)-\H^{B) - a^^dz. (4.5) 



Each resolvent has a polar integral kernel with a l/|x — x'| singularity, and 
the product of two resolvents will have a continuous kernel. In fact we can 
repeat this trick and obtain products of as many resolvents as we want, 
thus further improving the regularity of the integral kernel. Technical 
details will be given elsewhere. Actually this kind of argument can be 
used to show that all operators defined by integrals over complex contours 
have jointly continuous integral kernels. 

3. Although the order of limits in l|4.3(l is important for the argument below, 
it might be possible (at least under additional conditions on the spectrum 
of Hoo{B)) to interchange the order of limits. The important fact is that 
the thermodynamic limit has to be taken first: great care is to be taken 
when defining currents in the static limit for finite systems (for a discussion 
of this point in a related context see [21] . 

4. The result is valid for arbitrary magnetic field B and establishes the con- 
nection between the Hall conductivity and the Faraday effect. However, 
the quantum Hall effect requires high magnetic fields while the Faraday 
effect is usually considered at low magnetic fields. 

Proof. We start from the conductivity in the thermodynamic limit as given by 
Theorem 13.11 



ai2{B,T,uj) 



where 



■E 



fFD{z)T,{z,uj)dz 



(x,x) dx, (4.6) 



S(z, Lu) Pi{B){H^{B) - z)-^ P2{B){H^{B) -z- uj)-^ + (z^z-uj). (4.7) 

Since we made the assumption that {di,d2) C p{Hoo{B), then for < 
^"^''^ the integral over z on the contour T/s^;^ can be replaced with the integral 
on the contour F^^^ U F^^^ where (see also and (IT^ ): 



F^ 



X ±id : a < X < di 
d2 - di 



d2 - di 



[J{a + iy: -d < y < d} 



iy : —d < y < d 



(4.8) 
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and 



= \ x±id: x>d2- ^i^^l 



[j{d^~^^+'y- -d<y<dy (4.9) 
Accordingly, one can rewrite ai2{B,T,ijj) as 

'^-(^'^'-) = -^^i^E[ (4.10) 



2T:m?oj\Vt\ 

' ' s— 1 



S](z, w)(iz 

+ / Ufd{z) - l)S(z, tj)dz + / fFD{z)^{z, Lu)dz } (x, x) dx. 



Note that since the singularities of fpoiz) lie on ^^+iy, y £ (— cxo, oo), one can 
take ^ independent of /? i.e. one can take d = ^^^i^ in H2.29|l . At this point 
we take the limit /3 ^ oo. Since on we have \fFD{z)\ < 2exp[-/3(a;- ^li^)], 
and on Tl^ we have that Ifpoiz) — 1| < 2exp[— ~2;)], the last two terms 
in H4.10(l vanish in the zero temperature limit (full details about the control of 
various integral kernels will be given elsewhere). Hence we get: 

a.(i.,T = 0,.)^-^-^t/l/,^(^'-^ (-'-)^- (4-11) 

s—l \ ^ ) s.s 

An application of the Cauchy residue theorem shows that the two terms of 
S(z,u;) (see (|4.7|l ') will combine in the above integral and give 

2 1 2 

ai2(i?,r = 0,c.) = --^— ^ (4.12) 



/ / Pi{B){HooiB) - z + uj/2)-^P2{B){H^{B) - z - uj/2)- 
Jn Ur 



(x, x)(ix. 



\d2~di\ 



where F is any finite contour such that F C p{Hoo{B) + uj) for all \uj\ < 
and only enclosing the spectrum of Hao{B) below Now the integrand 

in (|4.12|l is analytic in in a neighborhood of the origin. By expanding the 
resolvents one obtains: 



ai2(B,r = 0,c^)= (4.13) 
-^F^E / (/ -Pi{B){H^{B)-z)-^P2{B){H^{B)^zr^ 

+ i [P^{B){H^{B) - z)-^P2{B){H^{B) - z)-2 

~Pi{B){HUB) ~ z)-^P2{B){H^{B) - z)-^] ] {^,^)d^ + 0{lo). 

J S.S 
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Apparently we have a first oder pole at a; = 0. But wc now prove that the 
singular term in the r.h.s. of (|4.13() is identically zero. Namely (when no spin 
variables appear the integral kernels below have to be understood as matrices 
in the spin space): 



1^ Pi{B)(H^{B) - z)-^P2{B){H^{B) ~z)-Hz^ (x,x) = 0. 



(4.14) 



Using the magnetic perturbation theory and the trick from H4.5|) one can prove 
that even though the integrand in H4.14|l has a quite singular kernel, after inte- 
gration with respect to z one gets a smooth kernel, exponentially localized near 
the diagonal (details will be given elsewhere). 

Let us notice an operator equality which makes sense on compactly sup- 
ported functions: 

-{H^{B) - z)~'P2{B)(H^{B) - z)-i = z[X2, {H^{B) - z)-% (4-15) 
m 

because the resolvent {Hoo{B) — z)^^ sends compactly supported functions into 
exponentially decaying functions (see (|3.7|l ). which are in the domain of X2- In 
fact, the operator on the right side has a nice integral kernel, given by 

{i[X2,(Hoo(S)-z)-i]}(x,y) =^(a;2-y2)Goo(x,y;z), (4.16) 

which is no longer singular at the diagonal and still exponentially localized near 
the diagonal, thus defining a bounded operator on the whole Hilbert space. 
After integration we get: 



2n 



J^Pi{B){Hoo{B) - z)-^P2(B){H^{B) - zY^dz = (4.17) 

'-^ J^P,{B)[{H^{B) - z)-\X2\dz - m[Pi(i3)n^,X2], 

where we used that Pi{B) and X2 commute. Note that the magnetic pertur- 
bation theory states that the integral kernel of Pi{B)Ii^ is smooth and expo- 
nentially localized near the diagonal. Therefore [Pi{B)n^,X2] wiU have the 
integral kernel: 

{[Pi(S)n^,X2]} (x,y) = (2/2 -:r2){Pi(S)n^}(x,y) (4.18) 

which is identically zero at the diagonal and proves 14.14|l . We can conclude at 
this point that: 



lim (712(5, T = 0,w) = (4.19) 

-=^E / {Pi{B){HUB)-z)-'P2{B){HUB)-z)-' 
I I Jn 

^Pi{B)iH^{B) - z)-^P2{B){Hoo{B) - z)-i | (x,x)dx. 

J s,s 

Now consider the r.h.s. of (|4.3|) . Since the magnetic field multiplying the 
spin will not change, our notation will only refer to i?2- Due to the stability 
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of the spectrum against small variations of the magnetic field, for sufficiently 
small AB, n^^+AS g^-j^ Q^i^ts and 



n 



B2+AB 



= ^ J^{iHoo{B2 + AB) - z)-i - (H^iB^) - z)-')dz. (4.20) 



By using the magnetic perturbation theory [231 with respect to AB (see also 
the discussion around 16.2|l ') one obtains: 

[n^^+^^-n^1(x,x) = ^ [ ( / dy{H^{B,) - z)-\^,y) (4.21) 
[Py(B2) • A(y - x)](i/oo(S2) - ^)"'(y,x)} dz + 0{{ABf). 

Now a very important identity is (see (|3.3() . 14.15|l and (14.161) ) 

A(y-x)(i7oo(B2)-z)-^(y,x) = (4.22) 
- A [{HUB2) ~ z)-^P(S2)(i?oo(S2) - 2)"'](y,x). 

The remainder in (AB)^ will have a smooth integral kernel after the integration 
with respect to z, hence we obtain: 

■ [ If \{H^{B)-z)-^Pi{B){H^{B)^z)-^P2{B){H^{B)-z)-^ 

- (i?oo(S) - z)-^P2{B){H^(B) - z)-^Pi{B){H^{B) - z)-i)dz] } (x,x). 

From l|4.23(l and H4.19|l we see that H4.3|l follows if we can prove that one can 
circularly permute the operators under the integral sign in (|4.23|) . One can prove 
this by interpreting 14.2311 as the thermodynamic limit of the corresponding 
expression on finite volume and then using the invariance of the trace under 
cyclic permutations. Alternatively one can prove it directly and in what follows 
we outline the proof. 

Due to the smoothing effect of the integral with respect to z, we can al- 
ways restrict ourselves to considering a product of only two integral opera- 
tors which commute with the discrete magnetic translations, and have kernels 
e^H.i^.y) K^{y^,y) and e*'"^''(y'^')i^2(y, x'). We therefore look at an absolutely 
convergent integral of the form (the anti-symmetric magnetic phases disappear 
when we look at the diagonal, see l|3.1(l ): 



/ dx / dyi^i(x,y)i^2(y,x) 

JO. JR3 



(4.24) 



with 

-ft^i,2(x,y) = i^i,2(x + 7,y + 7). 
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Then 



/ dx / dyXi(x,y)i^2(y,x) = V dx dyi^i(x, y + 7)i^2(y + 7, x) 

= V / dx / dyi4:i(x-7,y)ii:2(y,x-7) 
-g^Jo Jo 

= f dy f dx f i^2(y,x)ifi(x,y) (4.25) 

Jo JR3 Jr^ 

which gives the needed "trace cychcity" and the theorem is proved. □ 

We now turn to the question whether the hmit in 1)4.3(1 actually vanishes as 
is suggested by some heuristic arguments (see e.g. |SI])- We start by recalling 
some results about Wannier functions. Let ao{Bo) be an isolated part of the 
spectrum of Hoo{Bq) and fl^" the corresponding spectral projection. We say 
that n^" has a basis of exponentially localized (magnetic) Wannier functions 
if there exist a > 0, Wj G L2(]R^) © L2(R^), j = 1, 2, ...,p < oo satisfying (we 
denote by Wj{x.^ s) G C the values of Wj in x G M^^ and s G {1, 2}) 



2 p 

y/ |wj(x,s)pe2"l''ldx < M < oo, 



(4.26) 



such that the set of functions jwo ^} . -. „ with 
u;^- t;(x, s) = e'^'^-^^^T'^Wj- (x - 7, s) 

is a basis in the range of the projection {L'^ (M.^) © L^^jjs-j-j jj- ^j^g ^^^^ jg 
neglected, it has been proved in (22 that the existence of bases of exponen- 
tially localized Wannier functions is stable against small values of the magnetic 
field (i.e. Bq — 0). More precisely, if co is an isolated part of the spectrum of 
—/S. + V and the corresponding subspace has a basis of exponentially localized 
Wannier functions then, for sufficiently small B, a-o{B) is still isolated and the 
corresponding spectral subspace has a basis of exponentially localized magnetic 
Wannier functions. The methods in together with the magnetic perturba- 
tion theory [HI, HH, EH allow one to generalize the above result to arbitrary 
Bq and presence of the spin (as far as the spin-orbit term is sufficiently small) 
|12j . Now the existence of exponentially localized magnetic Wannier functions 
for an isolated part of the spectrum and for the value of the magnetic field B2 
in an interval around Bq allows one to write: 

J2 I n^..s(x,x)dx= [ y y |u;,-^(x,s)pdx 

= EE / rfxK(x-7,.s)p = f]^ / dxK(x,.)p=p. (4.27) 

. - „_i Jo 1 JR^ 



3,1 s=l •'^ j=l s=l ■ 



Thus the integrated density of states corresponding to the Fermi projection is 
constant in B2 in a small interval around Bq, hence this band gives no contri- 
bution in the r.h.s. of H4.3|l . 

For small fields, the above discussion can be summarized in: 
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Theorem 4.2. Suppose {di,d2) C /o(— A + V), d-z > di, and that the spec- 
tral subspace corresponding to {—oo,di] admits a basis of exponentially local- 
ized Wannier functions. Suppose that the spin-orbit interaction (see H2.2|l ) is 
small enough such that as (? decreases from oo to its actual value, we have that 
^^4^ S P(^foo(0)). Then for sufficiently small B: 

lim lim cri2(B. T, uj) = 0. (4.28) 



In particular all the derivatives of ai2{B,0) vanish for B = 0, and this 
substantiates Roth's result (formula (50) in [2]) for the first order correction 
in B at zero frequency. 

5 A closed formula for free electrons 

If y = it turns out that the conductivity tensor can be explicitly computed 
for all values of B and lu. The formula does not depend on whether we work in 
two or three dimensions. More precisely, we will show in this section that 

e'^n B 

ai2{B,uj)^— ^5^, (5.1) 

where n = n{T, /i, B) is the grand-canonical density. The formula H5.1|l is well 
known in classical physics and goes back at least to Drude but we are not 
aware of a known fully quantum derivation. The coincidence of classical and 
quantum formulas can be understood taking into account that the Hamiltonians 
involved (choose the symmetric gauge) are quadratic and it is known that for 
this class of operators classical and quantum computations coincide in many 
instances. While it is possible to derive H5.1|l by using the explicit form of 
the Green function or alternatively of eigenvalues and eigenprojections for the 
Landau Hamiltonian (see e.g. ^H]) we shall obtain it below only using resolvent 
and commutation identities. 

Let us only notice that when = we re-obtain formula (18) in 34 , while 
for a fixed frequency we get 



dB 

which is "the high frequency limit" or what Roth also calls "the free electron 
Faraday effect" in formula (51) from 31 . 

We begin by listing a few identities which are valid for a free electron on the 
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entire space. 



*[Pi(B),P2(i?)]-^ 

m c 

i[H^{B),P2iB)] = — Pi(P), 
m c 

32„2 



[ff^(B),[i/oc(i?),Pi(P)]] = ^PiiB), (5.2) 

[i/oo(P),[i?oo(i?),P2(P)]] = ^P2(P). 

Next, since in this case ^^^(x, x) does not depend upon x one has 

«(^'-) = -iE (5-3) 
1^ dzfFDiz) [Pi{B){HUB) - z)-^P2{B){H^{B) - z - co)-^ 
]}(0,s;0,s). 



Z —> Z — LU 



Commuting {Haa{B) — z — uo) ^ with P2(P) in the first term, and Pi{B) with 
{Hoo{B) — z + Lu)~^ in the second one, we obtain 

1^ d^/FcW [Pl(P)(i/oo(P)-2)-l(ifoc(P)-2-t^)-lP2(P) 

+ Pi{B){H^{B) - z)-\HUB) - z - Loy' 

■ [HUB),P2{B)]{H^{B) -z-uj)-^ 

+ (i?oo(P) - z + uj)-^Pi{B)P2{B){Hoo{B) - 

+ (P-OO(P) - 2 + C^)"'[-ffoo(P), Pl(P)] 

• (i/oo(P) - z + L^)-ip2(P)(i/oo(P) - z)"'] } (0,s;0,s) 
= I + II + III + IV. (5.4) 

Now I+III can easily be computed. Indeed, by cychc permutations (see (14. 25(1 ') 
one can cluster the two resolvents and then by the resolvent identity 

{A - z,)-\A - Z2)-^ = (zi - Z2)-^[{A - zi)-i - (A - Z2r% (5.5) 

one obtains four terms. Two of them vanish after the integration over z due to 
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the analyticity of the integrand while the other two give 

I + ///=— 5] (5.6) 

s=l 

1^ dzfFD{z)^[P2{B),Pi{B)]{H^{B) - z)-^Yq,s;6,s) 

2 

= — E{/f^^(^oo(i?))}(0,s;0,s) —n{T,fi,B). (5.7) 

Co' — 

s — 1 

In an analogous manner 

2 

III + IV =^Y. \ I '^'^^PD {z) (5.8) 



• {(iJoo(S) - z)-M^oo(S),P2(S)](iJoo(S) -Z-Uj)-^Pi{B) 

- (i7o,(B)-z)-i[i/oo(S),Pi(S)](i?oo(S)-z + u;)-iF2(S)} 
}(0,s;0,5). 

At this point we commute [iJoo(P), P2{B)] with (HooiB) — z — uj)~^ in the first 
term, and [Hoc{B), Pi{B)] with {Hoo{B) — z + Ld)~^ in the second term and use 
again. Some of the terms vanish after performing the integration over z 
and the remaining ones write as: 

- -{HUB) - z)-^[Hoo{B),P2{B)]Pi{B) (5.9) 

UJ 

~ -{HUB) - z)-^[HUB), [HUB),P2{B)]]{HUB) ~ z - lo)-^Pi{B) 

UJ 

- -{HUB) - z)-^[HUB),Pi{B)]P2{B) 

UJ 

- -{HUB) - z)-^[HUB), [HUB),Pi{B)]]{HUB) - z + lo)-^ P2{B). 

UJ 

Taking into account (|5.2I) the first and the third terms in H5.9II combine to 
--{HUB) - z)-^[HUB),Pi{B)P2{B)] 

which after integration over z is proportional to 

fFDiHUB))i[HUB),Pi{B)P2{B)] (5.10) 

= — {fFD{HUB))Pl{B)' - fFD{HUB))P2{B)^} , 

m c 

where we used the second and third identities in H5.2|l . Consider the uni- 
tary operator U which implements the coordinate change {U f){xi, X2, x^) — 
f{x2,-Xi,X3). Then one can prove that UPi{B)U* = -P2{B), UP2{B)U* = 
Pi{B) and UHUB)U* = HUB). This implies that 

UfFD{HUB))Pi{B)^U* = fFD{HUB))P2{B)^. 
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Since both operators have a smooth integral kernel, and because the rotation 
with U does not change the diagonal value of the integral kernel on the left hand 
side, it means that the contribution given by (|5.10(l is zero. 

Therefore we only remain with the second and fourth terms in H5.9|l . Using 
(|5.2|l . they become: 

- -^^{H^B) - z)-^P2{B){H^{B) - z - uj)-^Pi{B) 



- -K^{Hoo{B) - z)-^Pi{B){H^{B) - z + u)-^P2{B). (5.11) 

Using once more the cyclicity of the trace and comparing with the starting point 
1)5 .^ll . we obtain the remarkable identity 

B'^e^ 

II + IV= „ „ , a{B,uj). (5.12) 

Putting together l|5.4(l . 1)5. 6)1 . and ()5.12)l . we obtain the equation: 

B e S2g2 
a{B,uj) = n-\ a{B,uj), 

which gives ()5.1)l (see p.lO)! ^. 

6 Magnetic perturbation theory and the hnear 
term in B 

When 1/ it is no longer possible to obtain a closed formula for ai2{B,uj). 
Since in most physical applications the external magnetic field can be considered 
weak, an expansion in B up to the first or second order would be sufficient. In 
this section we show that aL{B, to) has an expansion in B to any order and write 
down the expressions of the first two terms. The first one gives the transverse 
conductivity at zero magnetic field and the second which is linear in B provides 
the Verdet constant. From 1)3.6)1 and ()3.8)) (in what follows by tr we mean the 
trace over the spin variable): 

aL{B,uj) = - ^ / dxltr / dz/poiz) (6.1) 
2vr|AL| Jaz, I Jr^ 

■ [ du{[(Px.i(0)-6Ai(x-u))XL(x,u;z)] 

■ [(-Pu,2(0) - 6A2(u - x'))i^L(u, x'; Z + LJ)] 
+ [(^x4(0) - &Ai(x - u))i^L(x, u; z - cj)] 

■ [{PuAO) - bA2{u - x'))i^L(u,x'; z)]}} 

Let us mention here that one cannot interchange the order of the above 
integrals. First one performs the integral with respect to u, then the integral in 
z, then we can put x = x' since the resulting kernel is smooth, and finally one 
integrates with respect to x over A^. 

When considering the expansion in boiaL{B,uj) we are left with the problem 
of the expansion of K\j^ (x, y; C) . This expansion is provided by the magnetic 
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perturbation theory as developed in [231 ■ Following the steps in ,23, in the case 
at hand one obtains: 



XL(x,y;z) G'^°^(x,y;z 



(6.2) 



-/ Gf(x,u;z) P„(0).A(u-y)Gr(u,y;z) 



(0), 



du 



+ f 4°) (x, u; z)a3Gf (u, y; z)du + 0{b') 

= Gf (X, y; .) + 6Gi°'^^") (x, y; z) + bd^-^ (x, y; z) + 0{h\ 

The above integrands are matrices in the spin variable, that is why the spin does 
not appear explicitly. The error term 0(6^) can also be fully controlled with the 
magnetic perturbation theory (actually arbitrary order terms can be computed; 
see 122] for details). Plugging the expansion (|6.2|l into H6.2|) and collecting the 
terms of zero and first order one obtains 



where the zeroth order term is: 

0^(0, cu) = - „ \ I / dxltr [ dzfpDiz) 

■ {Pi(0)(i/L(0) - z)-iP2(0)(i?L(0) - z + uj)-^ + {z^z- c^)}} 
while the first order correction reads as: 

/ \ orbit I \ 1 spin / \ 

aL,i(w) = a^^i (w) +0/^1 (w), 

where 

1 



(6.3) 
(6.4) 



orbit/ \ _ 



27r|A 



L\ Ja 



dx < tr 



dzfFoiz) 



(6.5) 
(6.6) 



du{- U(x-u)Gf (x,u;z)] [Pu,2(0)Gf (u,x';z + tj) 

P^Amf (x, u; z)] [A2{xi - x')Gf (u, x'; z + lo) 
P^Amf''^ (x, u; z)] [Pu,2(0)Gf (u, x'; z + 
(0)Gf (x, u; z)] [Fu,2(0)Gi°'"''") (u, x'; z + 
+ {z->z-uj)\\ 

du { [Px.i(0)Gi^-P™) (x, u; z)] [Pu,2(0)Gi°^ (u, x'; z + u) 

Px4 (0)Gf (x, u; z)] [Pu,2(0)Gi^^*") (u, x'; z + 
+ (z^z-t^)|| 

X — x' 



(6.7) 
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Now consider the expression A(x — y)G'^^ (x, y; z) appearing in the formula for 
«L,i('^)- Observing that it represents a commutator (see (|3.3() ) one has the 
identity 



A(x - y)Gi"^ (x, y; z) = ( A (x - y) ) G^^^ (x, y; z) 



1 



-naA [X,(i/L(0)-z)-i] ) (x,y) 



2m 



{(HLiO) - z)-\n3 A P){Hl{0) - z)-i}(x,y), 



(6.8) 



where X denotes the muhiphcation operator with x. By a straightforward (but 
somewhat tedious) computation one arrives at: 



/ \ orbit.l/ \ , orott'Z/ \ , spin/ \ 

o-lA^) = o-L.i ' (^)+aL,i (w) + aL,i (w) 



bit, 2/ 



where 



orbit, 1 ( \ 
Oil = 



dx < tr / dzfpDiz) 



(6.9) 



(6.10) 



^ P„(0)(i/L(0) - z)-ip„(0)(ffL(0) - z - u:)-^ 

a = l 
2 

J2 PAOKHl{0) - z)-^PamHL{0) - z + c^)-i 
=1 

2 

J2 Po.mHL{0) - z)-^P^{0){Hl{0) - z)~ 



2 



x,x), 
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W(^^l(0)- 


z)-^Pi(0)(Hl(0) - 


2)"' 

/ 


■ ^2(0)( 


-'-'V / 




^P2(0)(Hl(0) - z- 


/ 


+ Pi{0 


){Hl{0) 




-^P2(0)(Hl(0)- z) 


-1 


•Pi(0) 


Hl(0) - 




^P2(0)(HLiO) - z- 


/ 


- Pi{0 


)(^l(0) 


— z - 


f w)"lPl(0)(i?L(0) 




■ P2i0) 


Hl(0) - 

-'--'V / 


- z + 


ujy'^P2io)(HLio) - 


/ 


+ Pi{0 


)(^l(0) 


— z - 


f a;)"ip,(0)(i?L(0) 




•Pi(0) 


Hl{0) - 


- z + 


tj)"^P2(0)(i/L(0) - 


/ 


-Pi(0 


)(^l(0) 


~z) 


-^P2{0){Hl{0)-z 




•Pi(0) 


^^l(0) - 


- z — 


C^)-lP2(0)(i/L(0)- 


-z-Luy^ 


+ Pi(0 


){HLiO) 


^z) 


"'P2(0){7?l(0)-z 




■P2(0)( 


Hl{0) - 


- z — 


c^)-1Pi(0)(Hl(0)- 


-z-ujy^ 


-Pi(0 


){Hl{0) 


— z - 


f c^)-ip2(0)(i/L(0) 




•Pi(0) 


HlIO) - 


-z)- 


iP2(0)(i/L(0)-z)- 


1 


+ Pi(0 


)(^l(0) 


— z - 


f L^)-lP2(0)(i/i(0) 


-z)-' 


•P2(0) 


{Hl{0) 


-z) 


-iPi(0)(i/,.(0)-z) 


-'}}(x,x) 



and 

<r(^) = -^f^ ^^{tr / (6.12) 

• {[Pl(0)(i/L(0) - z)-V3(i?L(0) - z)-lP2(0)(i/L(0) - 2 - C.)-!] 

+ [Pi(0)(i7L(0) - z)-ip2(0)(i/L(0) - z - c^)-V3(i/L(0) - z - u^r^] 
+ [Pi(0)(i/L(0) - z + Luy'asiHLiO) - z + L^)-ip2(0)(i/i(0) - z)-'] 

+ [Pi(0)(i7L(0) - z + c^)-ip2(0)(i/L(0) - z)-V3(i?L(0) - z)-i] }} (x,x). 

7 The periodic case 

Now consider the case when V is periodic. In this case, after taking the ther- 
modynamic hmit one can replace (see (IS-BI) ) J^^ with J^^ and rewrite 
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as integrals over the Brillouin zone 

2 

• + ^a)(/i(k) - zy^ip^ + ka){h{k) -z- uj)-^ 

2 

+ + ^a)(/j(k) - ^)"^(Pa + fcc)(/l(k) - Z + 

Q=l 

2 

- + fc„)(/i(k) - z)-ip„(0)(/i(k) - z)-i } (x,x), 

^^'^ - ^\ L L 1 '^^^^^^^ ^'-'^ 

• {-(pi + fci)(/i(k) - + fci)(;i(k) - 



•(P2 + 


fc2) 


\m - 




-\V2^k2){h{k)-Z- 


-)-^ 




+ /ci 


)(Mk)- 


-A 


"'(p2+fc2)(/l(k) -Z) 


-1 




fcl) 


:m - 




"'(p2+fc2)(Mk) - 




-(Pl 


+ fci 


)(Mk)- 


- z 


+ c^)-i(pi + fci)(Mk) 


- Z + LJ)-^ 


•^2 + 


fc2) 


:/j(k) - 


z 4- 






+ (Pl 


+ fci 


)(Mk)- 


- z 


+ L^)-l(p2+fc2)(Mk) 


- z+ujy^ 


•(Pi4 


fcl) 


;/i(k) - 


z 4- 


-^)"'(p2+fc2)(Mk) - 




-(Pi 


+ fcl 


)(Mk)- 


-A 


"'(p2+fc2)(Mk)-2 




•(Pi4 


fcl) 


;/i(k) - 


z - 


-^)-\p2 + h2){h{V) - 


-z-oj)-^ 


+ (Pi 


+ fcl 




-A 


-\p2+k2){h{k)-Z 




•(P2 + 


fc2) 


\m - 


z - 


- uj)-\p, + h){h(k) - 


-Z-LJ)-^ 


-(Pl 


+ fcl 


)(Mk)- 


- z 


+ Uj)-\p2+k2)(h(k) 




•(P14 


fcl) 


\m - 




-\p2 + k2){h{k)~z)- 


1 


+ (P1 


+ fcl 


Mm - 


- z 


+ Uj)-\p2+k2){h{k) 




• (P2 + ^2) 


:/i(k) - 




-i(pi + fci)(/i(k) -z)- 


}(x,x) 



<i (^) £ i {tr dzfM^) (7.3) 

• { [(pi + fcl)(/l(k) - Z)-V3(/I(k) - Z)-\P2 + k2){h{k) - Z - C^)-l] 

+ [{pi + fci)(/i(k) - z)-\p2 + k2){h{k) - z - uj)-'a3{h{k) - z - ^)-i] 

+ [{pi + k,){h{k) - Z + Uj)-^<J3{h{k) - Z + Uj)-\p2 + k2){h{k) - Z)-1] 

+ [{pi + fcl)(/l(k) - Z + U)-\P2 + k2){h{k) - Z)- V3(/i(k) - Z)-1] } } (X, X). 

Finally, for the convenience of the reader only interested in applying the 
theory to the case when one assumes that the Bloch bands and functions are 
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known (as for example from Kohn-Luttinger type models), we write H7.1|) - (|7.3|l 
in terms of Bloch functions and energies. The important thing here is that 
no derivatives with respect to the quasi-momentum appear. With the usual 
notation (here (, ) denotes the scalar product over the spin variables): 

7ry(a,k)= {ui{x,k),{pa +ka)uj{x,k))dx, (7.4) 

and after some rearrangements, the terms coming from the orbital magnetism 
are: 



orbit,!/ \ 



T^ilf rfkJ^fe(a,k)p/;^^(A,(k)) (7.5) 



,,2Y^I^ ^„ 1.^1 2 /FD(Aj(k)) - /i.Z5(Ai(k)) 



[(A,(k)-AKk))2-c.2](A,-(k)-A,(k)) 



= 2;;^W L ^ i /r '^"^^''^"^ ^^-^^ 

ni ,712 ,^3 ,n4>l ^ 
J "^714711 (1 1 ^)'^nin2 (^i k)'7r,-2,2n3 (2, k)7rn3r),4 (2, k) 

1 (z - A„,(k))(z - A„,(k))(z - X„,{k))iz+Lu- A„,(k)) 



'''ri4rii (Ij k)7r„j„2 (2, k)7r„2„3 (1, k^TTn^m (2, k) 



+ 



+ 



+ 



(z 


^ A„j 


(k))(^-A„, 
T^uiUi (Ij k) 


(k))(z-A„3(k))(z + t^-A„,(k)) 

^ni n2 (I7 k)7r^2n3 (2, k)7r7i37i4 (2, k) 






(z 


— to - 


An,(k))(z- 
^rii-ni (Ij k) 


- - A„2 (k))(z - - A„3 (k))(z - 

^"in2 (2j k)7T„2„3 (1, k)7T„3„^ (2, k) 


An4 


(k)) 


(z 


~ LU - 


A„,(k))(z- 

^riin-i (Ij k) 


- LO - A„2(k))(z - a; - A„3(k))(2 - 

^nin2 (2, k)'7r^2n3 (1; k)'7r^3^^ (2, k) 


An4 


(k)) 


(z 


— A„j 


(k))(z + ^- 

'''n4ni (Ij k) 


- A„2(k))(z +UJ - A„3(k))(z +UJ - 

^nin2 (2, k)'7r^2n3 (2, k)'7r^3^^ (1, k) 


An4 


(k)) 


(z 


— A„j 


(k))(z + ^- 


' A„2(k))(z +UJ - A„3(k))(z + w - 


An4 


(k)) 




7''n4rii (Ij k)7T„j„2 (2, k)7T„2Tl3 (1, k)7T„3„^ (2, k) 






(z 


— to - 


A„,(k))(z- 


^A„2(k))(z-A„3(k))(z-A„,(k)) 








'f^niiii (1 ; k)7r,ij„2 (2, k)';T,i2,j3 (2, k)7r,i3„^ (1, k) 






(z 


— — 


A„,(k))(z- 


A„2(k))(z-A„3(k))(z-A„,(k))^ 







As for the spin contribution, with the notation 

Sy(k) := / (Mi(x,k),o-3Wj(x,k))dx, (7.7) 
Jn 
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one has: 




{ 



(An2 (k) - z) (An3 (k) - z) ( (k) - z - uj) 



^711712 k)s722n3 (k) TT^^^ 7^ (2, k) 



+ 



(A„2 (k) - z) ( A„3 (k) - z - w) (A„^ (k) - z - w) 



+ 



(k) (2,k) 



(A„2 (k) - z + c^)(A„3 (k) - z + u){Xni (k) - z) 



+ 



^nin2 (l: k)7r7T,2n3 (2, k)^^^^^-!^ (k) 



(A„, (k) - z + w) (A„3 (k) - z) (A„, (k) - z) 



8 Conclusions 

We presented in the present paper a method which shed new hght on the quan- 
tum dynamics/optical response in bulk media in the presence of a constant 
magnetic field. We applied the gauge invariant magnetic perturbation theory 
and gave a clear and very general way of dealing with long range magnetic 
perturbations. 

The formal connection with the integer Quantum Hall effect was established 
in l|4.2|l . Equations (|6.9(l - (|6.12|l and (|7.4|I - H7.8|I contain our main result concern- 
ing the Verdet constant and the Faraday effect: it gives the linear term in B of 
the transverse conductivity in terms of the zero magnetic field Green function. 
They open the way of using the recently developed Green function techniques 
for the calculation of optical and magneto-optical properties of solids, to the 
case when an external magnetic field is present. Our method can be applied to 
ordered, as well as to random systems (with the appropriate average over con- 
figurations). Of course, in the last case one has to assume ergodicity properties 
in order to insure convergence of results in the thermodynamic limit. Layers or 
other geometries can also be considered. 

There are many subtle and difficult mathematical questions left aside in this 
paper, as those related to the thermodynamic limit, the convergence of infinite 
series over Bloch bands, and the low frequency limit when the Fermi energy lies 
in the spectrum. Another open problem is to consider self-interacting electrons 
and to investigate the exciton influence on the Faraday effect. These questions 
will be addressed elsewhere. 

Our results are not only theoretical. In a future publication we will use the 
residue theorem in equations H7.4|l - 17.8|l to calculate the Verdet constant for 
various finite band models, and compare our results with the existing experi- 
mental data. Moreover, our results will be shown to imply those of Roth |3J 
and Nedoluha HO]. 

9 Acknowledgments 

This paper has been written at the Department of Mathematical Sciences, Aal- 
borg University. G.N. kindly thanks the Department for invitation and financial 
support. 



26 



H.C. acknowledges support from the Statens Naturvidenskabelige Forskn- 

ingsrad grant Topics in Rigorous Mathematical Physics^ and partial support 
through the European Union's IHP network Analysis & Quantum HPRN-CT- 
2002-00277. 



References 

[I] Angelescu, N., Bundaru, M., and Nenciu, G.: "On the Landau diamag- 
netism", Commun. Math. Phys. 42, 9 (1975). 

[2] Banhart, J.: "Optical Conductivity of Disordered Alloys Calculated from 
First Principles", Phys . Rev. Lett. 82, 2139 (1999). 

[3] Bellissard, .J., and Nakamura, S.: "Low energy bands do not contribute to 
quantum Hall effect", Commun. Math. Phys. 131, 283 (1990). 

[4] Bennett, H.S., and Stern, E.A.: "Faraday effect in solids", Phys. Rev. 137, 
A448 (1965). 

[5] Boswarva, LM., Howard, R.E., and Lidiard, A.B.: "Faraday effect in semi- 
conductors" , Proc. Roy. Soc. (London) Series A. Mathematical and Physical 
Sciences, 269, 125 (1962). 

[6] Boswarva, I.M., and Lidiard, A.B.: "Faraday effect in semiconductors. H" , 
Proc. Roy. Soc. (London) Series A. Mathematical and Physical Sciences 278, 
588-609 (1964) 

[7] Bratelli, O., and Robinson, D.W. Operator Algebras and Quantum Statistical 
Mechanics 2 (Springer- Verlag, Berlin, Heidelberg, New York, Second Edition 

1997) 

[8] Bouclet, J-M, Germinet, F., Klein, A., and Schenker, J.H.: "Linear response 
theory for magnetic Schrdinger operators in disordered media", J. Funct. 
Anal. 226 no.2, 301 (2005). 

[9] Butler, W.H.: "Theory of electronic transport in random alloys: Korringa- 
Kohn-Rostoker coherent-potential approximation", Phys . Rev. B 31, 3260 
(1985). 

[10] Combes J.M., and Thomas, L.: "Asymptotic behaviour of eigenf unctions 
for multiparticle Schrodinger operators", Commun. Math. Phys. 34, 251 
(1973). 

[II] Cornean, H.D.: "On the magnetization of a charged Bose gas in the canon- 
ical ensemble", Commun. Math. Phys. 202, 1 (2000). 

[12] Cornean, H.D., Nenciu, A., and Nenciu, G., in preparation. 

[13] Cornean, H.D., and Nenciu, G.: "On eigenfunction decay for two- 
dimensional magnetic Schrodinger operators", Commun. Math. Phys. 198, 
671 (1998). 

[14] Cornean, H.D., and Nenciu, G., in preparation. 



27 



15] A. Gonis, Green Functions for Ordered and Disordered Systems (North- 
Holland, Amsterdam, 1992). 

16] Huhnc. T., and Ebcrt, H.: "Fully relativistic description of the magneto- 
optical properties of arbitrary layered systems", Phys . Rev. B 60,12982 
(1999). 

17] Jimenez-Gonzales, H.J., Aggarwal, R.L., and Favrot, G.: "Infrared Faraday 
rotation of n-type InSb", Phys. Rev. B 49, 4571 (1994). 

18] Joynt, R., and Prange, R.: "Conditions for the quantum Hall effect", Phys. 
Rev. B84, 3303 (1984). 

19] Luttingcr, J.M.: "The Effect of a Magnetic Field on Electrons in a Periodic 
Potential", Phys. Rev. 84, 814 (1951). 

20] Nedoluha, A.: "Perturbation-Theoretical Approach to Magneto-Optical 
Phenomena in the Spectral Region of Zero Absorption", Phys. Rev. 139, 
A1489 (1965). 

21] Nenciu, A., and Nenciu, G. : "Existence of exponentially localised Wannier 
functions for nonperiodic systems", Phys. Rev. B 47, 10112-10115 (1993). 

22] Nenciu, G.: "Dynamics of band electrons in electric and magnetic fields: 
rigorous justification of the effective Hamiltonians" , Rev. Mod. Phys. 63, 91 
(1991). 

23] Nenciu, G.: "On asymptotic perturbation theory for quantum mechanics: 
Almost invariant subspaces and gauge invariant magnetic perturbation the- 
ory", J. Math. Phys. 43, 1273 (2002). 

24] Nockcl, J.U., Stone, A.D., and Baranger, H.U.: "Adiabatic turn-on and 
the asymptotic limit in linear-response theory for open systems" , Phys. Rev. 
B 48, 17569-17572 (1993). 

25] Pedersen, T.G.: "Tight-binding theory of Faraday rotation in graphite", 
Phys. Rev. B 68, 245104 (2003). 

26] Pcicrls, R.E.: "Zur Theorie des Diamagnetismus von Leitungelectronen" , 
Z. Physik 80, 763 (1933). 

27] Perlov, A., Chadov, S., and Ebert, H.: "Green function approach for the 
ab initio calculation of the optical and magneto-optical properties of solids: 
Accounting for dynamical many-body effects", Phys . Rev. B 68, 245112 
(2003). 

28] Pidgeon, C.R., and Brown, R.N.: "Interband Magneto-Absorption and 
Faraday Rotation in InSb", Phys . Rev. 146, 575 (1966). 

29] Roth A. P., and Fortin, E.: "Magneto-optics in semiconductor crystals of 
arbitrary orientation" , Phys. Rev. B 18, 4200 (1978). 

30] Roth, L.M. : " Theory of Bloch electrons in a magnetic field" , J. Phys. Chem. 
Solids 23, 433 (1962). 



28 



[31] Roth, L.M.: "Theory of Faraday effect in sohds", Phys. Rev. 133, A542 
(1964). 

[32] Schwingcr, .J.: "On Gauge Invariance and Vacuum Polarization", Phys. 
Rev. 82, 664 (1951). 

[33] Sondheimer, E.H., and Wilson, A.H.: "The theory of the magneto- 
resistance effects in metals", Proc. Roy. Soc. (London) A 190, 435 (1947). 

[34] Streda, P.: "Theory of quantized Hall conductivity in 2 dimensions", J. 
Phys. C.-Solid State Phys. 15, L717 (1982). 

[35] Thouless, D.J.: "Wannier functions for magnetic sub-bands", J. Phys. C.- 
Solid State Phys. 17, L325 (1984). 



29 



